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SUMMARY 


OBJECTIVE 

Determine  the  feasibility  of  using  an  optical  system  composed  of  a  remote  fiber-optic  probe  coupled  to 
a  scanning  partial  coherence  Michelson  interferometer  for  localized  thickness  measurement  of  thin  silicon 
films. 


RESULTS 

A  partial  coherence  interferometer,  combined  with  an  optical  fiber  coupler  probe  has  been  designed, 
assembled,  characterized,  and  used  to  accurately  measure  the  thickness  of  a  thin  silicon  substrate 
charge-coupled  device  (CCD).  The  system  has  been  configured  in  a  way  that  it  can  be  adapted  to  measure 
thick  samples,  and  to  perform  in  situ  measurements  of  samples  during  the  etching  process.  Feasibility  of 
the  measurement  technique  has  been  demonstrated  on  actual  processed  CCD  devices  with  a  nominal 
thickness  of  8  pm,  and  an  estimated  precision  of  1.3  percent,  or  0.1  pm,  has  been  achieved. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A  remote  fiber-optic  probe  coupled  to  a  scanning  Michelson  interferometer  can  be  configured  into  an 
optical  system  to  perform  in  situ  precision  measurements  of  thin  silicon  film  thickness. 

Such  a  system  can  function  in  realtime  during  a  manufacturing  process,  such  as  film  deposition  or 
removal,  and  can  be  used  to  aid  in  the  control  of  the  process. 

The  prototype  laboratory  model  developed  in  this  program  should  be  modified  to  demonstrate  this 
capability  and  assess  its  limitations. 
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INTRODUCTION 


For  a  particular  Navy  application,  it  is  necessary  to  “thin”  the  substrate  of  a  silicon  charge- 
coupled  device  (CCD)  after  the  device  sensor  array  and  microcircuitry  have  been  processed.  This 
enables  back-plane  illumination  of  the  CCD  without  significant  attenuation  or  scattering  of  the  optical 
signals  incident  on  the  sensor.  Thinning  is  done  by  means  of  an  etching  technique,  which  must  be 
controlled  precisely,  and  therefore,  it  is  advantageous  to  be  able  to  measure  the  substrate  thickness 
during  the  etching  process.  Furthermore,  the  final  thickness  of  the  etched  silicon  substrate  must  be 
known  precisely,  with  an  accuracy  on  the  order  of  0.1  pm. 

A  partial  coherence  interferometer,  combined  with  an  optical  fiber  coupler  probe  has  been 
designed,  assembled,  characterized,  and  used  to  accurately  measure  the  thickness  of  a  thin  silicon 
substrate  CCD.  The  system  has  been  configured  in  a  way  that  it  can  be  adapted  to  measure  thick 
samples,  and  to  perform  in  situ  measurements  of  samples  during  the  etching  process.  Feasibility  of  the 
measurement  technique  has  been  demonstrated  on  actual  processed  CCD  devices  with  a  nominal 
thickness  of  8  pm,  and  an  estimated  precision  of  1.3  percent,  or  0.1  pm  has  been  achieved. 


PARTIAL  COHERENCE  INTERFEROMETRY 

Partial  coherence  exists  in  an  interferometric  system  whenever  the  optical  source  cannot  be  con¬ 
sidered  purely  monochromatic.  That  is,  it  has  a  spectral  distribution  with  a  bandwidth  (in  wavenum¬ 
bers)  comparable  to  the  inverse  of  the  differential  path  length  of  the  optical  beams  propagating  in  the 
interferometer.  The  spectral  distribution  is  characterized  by  a  central  wavelength,  A0,  and  a  spectral 
bandwidth,  <5A,  or  alternatively,  in  wavenumbers  by  a0  =  \/kQ  and  da  =  6X/Xl  (see  figure  1).  The 
coherence  length,  lc,  of  the  source  is  a  measure  of  the  distance  over  which  two  beams  or  wave  pack¬ 
ets,  which  are  derived  from  the  source  by  means  of  amplitude  division,  will  produce  high-visibility 
interference  fringes  when  they  are  recombined.  The  expression  for  lc  in  terms  of  the  spectral  parame¬ 
ters  is' 


lc=  \/do  =  Xl/dX  (1) 

The  relationship  between  source  bandwidth,  coherence,  and  fringe  visibility  is  shown  schematically 
in  figure  2  where  the  output  from  the  source  is  represented  by  a  wave  packet  of  length  lc.  The  wave 
packet  is  divided  into  two  equal  amplitude  waves,  then  recombined  after  traversing  optical-path  lengths 
differing  by  an  amount,  A p  •  (An  optical  system  that  can  perform  this  function  is  a  Michelson  inter¬ 
ferometer  described  later.)  In  figure  2a,  the  waves  are  assumed  to  have  traveled  over  optical  paths 
differing  by  an  amount,  A p,  which  is  greater  than  the  coherence  length  of  the  source.  In  this  case, 
there  is  no  longitudinal  spatial  overlap  of  the  wave  packets  when  they  are  recombined,  and  the  beams 
are  incoherent  with  an  intensity  equal  to  the  sum  of  the  individual  wave  intensities.  If  the  waves 
traverse  optical  paths  differing  by  an  amount,  A p  <lc>  then  there  is  longitudinal  spatial  overlap  and 
the  waves  are  partially  coherent  when  they  are  recombined.  (Complete  coherence  exists  only  for 
A p  =  0.)  In  this  case,  shown  in  figure  2b,  the  intensity  depends  on  the  phase  difference  between  the 
waves  being  proportional  to  the  optical  path  difference. 

As  the  path  difference,  A p,  is  increased  from  0,  the  intensity  will  vary  in  a  cosinusoidal  way,  with 
a  continuously  decreasing  amplitude.  This  amplitude  factor  is  a  function  of  lc  and  also  of  A p,  and  is 
known  as  the  fringe  visibility2  or  mutual  coherence  function  7  3  This  function  is  determined  by  the 
spectral  distribution  of  the  source,  and  therefore,  source  spectral  characteristics  are  a  critical  factor  in 
the  ability  to  measure  sample  thickness  using  Michelson  interferometry. 
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SPECTRAL  DISTRIBUTION  WAVE  PACKET 


lc  =  l /6a  =  XI/6X 


Figure  1.  Optical  source  spectral  energy  distribution  and  wave  packet 
emission  representation  of  source  output.  Coherence  length  lc  of 
source  is  determined  by  spectral  width  Al¬ 


la)  A p>lc  I(Ap)  =  /1+/2 


(b)  A p<lc  I(bp)  =  +  +y(Ap)cos(2jto0Ap)] 

Figure  2.  Recombination  of  amplitude-divided  wave  packet  from  source, 
(a)  Path  difference  between  waves  A p>lc-  (b)  Path  difference  Ap  S /c  . 
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As  an  example  of  the  interference  of  waves  produced  by  a  source  with  an  extended  spectral 
distribution,  the  case  of  a  Gaussian  distribution  is  considered  here.  If  the  source  optical  field  is  repre¬ 
sented  by  the  function  Es(a),  then  the  field  produced  by  the  interfering  beams  for  a  wavenumber 
component  a  at  a  time,  t,  is 

E(o)  =  Re{jEs(o)(exp[i(2:topi  -  wt)}  +  exp[/(2;rqp2  -  a>t)])}  (2) 


where  P\  and  pz  are  the  optical  paths  of  the  two  waves.  The  components  of  the  field  are  assumed  to 
have  the  same  polarization  and  that  the  optical  paths  are  in  a  nondispersive  medium.  (The  more  gen¬ 
eral  case  including  a  dispersive  optical  path  is  discussed  in  r  ppendix  A.)  Also  assumed  are  the  source 
can  be  treated  as  a  point  source,  and  that  there  is  no  angular  misalignment  between  the  interfering 
waves.  Under  these  conditions,  there  is  no  transverse  spatial  variation  of  E(o)-  The  time-averaged 
intensity  for  a  particular  wavenumber  component  of  the  source  is  then  given  by 

J(o)  =  \E(o)\2  =  i|£s((7)|2[l  +  cos(2zroA/>)]  (3) 


where  A p=p\-pz-  The  total  intensity  is  found  by  integrating  the  intensities  resulting  from  each  com¬ 
ponent  of  the  source  spectrum.  The  total  time-averaged  intensity  is  then  given  by 

I,{Ap)  =  j  I  (a)  do  =  i  j  \Es(o)\2do  +  ij  \Es(o)\2  cos(2zioAp)do  .  (4) 

The  first  term  above  is  constant  and  represents  the  background  intensity  level  equal  to  half  the 
source  intensity,  or  equivalently,  /,(»).  The  second  term  varies  with  respect  to  the  path  difference, 

A p,  and  for  nondispersive  media  it  is  simply  the  Fourier  cosine  transform  of  the  source  energy  spec¬ 
tral  distribution. 

For  a  point  source  with  a  normalized  Gaussian  spectral  distribution  given  by 

|Es(eOI2  =  (l/<5a)  exp  {-  n[(,o  -  o0)  /6o]2}  (5) 


the  normalized  intensity  as  a  function  of  A p  is  from  equations  (4)  and  (5) 

It(Ap)  =  I  {1  +  exp[-zr(Apdo)2]  cos(2jro0Ap)}  (6) 


From  this  expression,  the  general  observations  made  regarding  the  interference  intensity  are 
confirmed.  First,  in  regions  of  nonzero  fringe  visibility  (i.e.,  A p<  l/<5a ),  the  intensity  varies  as  the 
cosine  of  a  term  proportional  to  the  optical-path  difference,  A p,  and  with  a  period  equal  to  the  mean 
wavelength  of  the  source  distribution,  X0  =  l/o0  ■  Second,  the  fringe  visibility  decreases  as  the  differ¬ 
ence  in  optical-path  length,  A p,  increases.  Also,  the  extent  of  the  fringe  visibility  is  correlated  to  the 
quantity,  lc,  the  coherence  length  of  the  source  equal  to  the  inverse  of  the  spectral  width  of  the 
source  lc  =*  1  /do  . 

Figure  3  is  a  plot  of  the  variable  intensity  term  in  equation  (6)  as  a  function  of  A p-  This  type  of 
graph  known  as  the  source  interferogram  can  be  produced  by  a  scanning  Michelson  interferometer. 
The  source  parameters  used  in  equation  6  to  generate  this  plot  are  k0  =  830  nm  and  dX  =  60  nm. 
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Figure  3.  Source  interferogram  calculated  for  a  Gaussian  spectral  distribution 
source.  Assumed  parameters  are  A0  =  830  nm  and  (5A  =  60  nm. 


Although  the  above  statements  regarding  the  features  of  the  interferogram  are  made  for  a  source 
with  a  Gaussian  distribution,  they  are  true  qualitatively  for  a  general  source.  In  practice,  the  source  is 
not  a  true  Gaussian,  in  that,  it  can  exhibit  structure  and  asymmetry  in  the  spectrum.  Asymmetry  in 
the  spectrum  has  the  effect  of  extending  the  width  of  the  visibility  function,  while  periodic  structure 
causes  sidelobes  of  appreciable  visibility.  These  effects  are  shown  in  the  plots  in  figures  4  and  5. 
Three  different  source-spectral-distribution  curves  are  shown  in  figure  4:  (a)  corresponding  to  a  pure 
Gaussian,  (b)  Gaussian  with  symmetric  periodic  structure  superimposed,  and  (c)  Gaussian  with  asym¬ 
metric  periodic  structure.  In  both  b  and  c,  the  structure  period  is  AA  =  10  nm,  or  in  wavenumbers, 
Act  =  0.0145  Pm"1  .  Each  distribution  has  the  same  total  spectral  width  of  <5A  =  40  nm  and  a  central 
wavelength  of  A0  =  830  nm.  The  corresponding  source  interferograms  for  positive  values  of  the  path 
difference  are  shown  in  figure  5.  As  seen  in  5b,  the  visibility  for  the  symmetric  distribution  with 
periodic  structure  is  significantly  increased  for  values  of  the  path  length  difference  in  the  region  of 
1/Act.  This  is  expected  since  the  interferogram  is  just  the  Fourier  cosine  transform  of  the  spectral 
distribution.  Similar  results  are  evident  for  the  asymmetric  case  in  figure  5c.  In  addition,  observe  that 
the  envelopes  near  0  and  1/Act  are  broadened  due  to  the  asymmetry.  Note  that  in  both  cases,  al¬ 
though  the  visibility  function  is  modified,  the  period  of  the  fringes  remains  at  the  mean  wavelength. 
These  effects  of  spectral  structure  and  asymmetry  must  be  taken  into  account  when  measurements  of 
film  thickness  are  made  using  partial  coherence  interferometry.  This  point  is  further  explored  in  the 
discussion  section. 
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Figure  4.  Spectral  distributions  for  (a)  a  pure  Gaussian  source, 
(b)  a  Gaussian  with  periodic  structure  superimposed,  and  (c) 
a  Gaussian  with  asymmetric  periodic  structure  superimposed. 


y 


l/A a  =  69  pm  A p  (pm) 


Figure  5.  Source  interferograms  calculated  for  the  source 
spectral  distributions  in  figure  4. 
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EXPERIMENTAL  CONFIGURATION 


The  properties  of  an  interferogram,  obtained  from  a  Michelson  interferometer  using  a  partial 
coherence  source,  can  be  used  to  measure  the  thickness  of  a  thin  silicon  CCD  sample.  Consider  the 
schematic  diagram  of  the  experimental  configuration  in  figure  6.  The  source  is  a  commercially  avail¬ 
able  edge-emitting,  light-emitting  diode  (ELED),  packaged  in  a  hermetic  housing,  and  coupled  to  a 
multimode  fiber  pigtail.  The  housing  contains  a  thermoelectric  cooler  and  a  thermistor  to  stabilize  the 
temperature  and  wavelength.  The  ELED  emits  in  the  wavelength  region  near  0.83  (im,  with  a  spectral 
width  of  approximately  40  nm.  The  exact  spectral  energy  distribution  depends  on  temperature  and 
operating  current.  Source  distribution  effects  will  be  discussed  in  a  subsequent  section. 


Figure  6.  Experimental  configuration  for  the  measurement  of  CCD 
sample  film  thickness  using  an  optical  fiber  coupler  probe,  a  Michelson 
interferometer,  and  a  partially  coherent  source. 


The  ELED  fiber  pigtail  is  spliced  to  one  input  port  of  a  multimode-fused  fiber  bidirectional  cou¬ 
pler.  The  coupler  is  a  3-dB  (or  50/50)  optical  power  splitter  at  0.83  nm.  One  of  the  coupler  output 
ports  is  positioned  near  a  graded  index  (GRIN)  rod  lens  where  its  output  can  be  either  collimated  or 
focused  onto  the  CCD  sample  to  be  evaluated.  The  second  fiber  output  can  serve  as  an  input  port  for 
a  visible  wavelength  (HeNe)  laser  used  to  align  the  interferometer  portion  of  the  measurement  system, 
or  it  can  be  coupled  to  a  spectrometer  for  monitoring  the  source  spectral  properties  during  a 
measurement. 


6 


The  optical  beam  incident  on  a  CCD  sample  of  thickness,  t,  is  partially  reflected  and  transmitted 
at  the  first  surface.  The  transmitted  beam  is  again  partially  reflected  at  the  back  surface.  The  front- 
and  back-surface-reflected  beams  are  coupled  back  into  the  GRIN  lens  and  the  fiber  coupler.  Use  of 
a  multimode  fiber  maximizes  the  input/output  coupling  efficiency  of  the  optical  beams.4’5  The  fourth 
coupler  port  transmits  the  sample-reflected  beams  to  the  Michelson  interferometer.  The  two  waves  are 
divided  by  the  Michelson  beamsplitter,  and  one  set  of  wave  packets  is  delayed  with  respect  to  the 
other  by  the  variable  path  length  difference,  2pt,  between  the  two  arms  of  the  interferometer.  This 
path-length  difference  is  controlled  by,  and  is  proportional  to,  the  voltage,  Vp,  applied  to  the  linear 
piezoelectric  transducer.  The  beams  are  recombined  at  the  beamsplitter  and  interfere  on  the  detector 
whose  output  current  is  proportional  to  the  intensity  of  the  interfering  beams.  This  detector  current  is 
amplified  and  the  amplifier  voltage  is  either  displayed  on  an  oscilloscope  as  a  function  of  Vp  (i.e.,  p<), 
or  it  can  be  converted  to  a  digital  signal  and  stored  in  a  computer.  A  linear-ramp  voltage  is  applied  to 
the  piezoelectric,  and  the  fringe  intensity  versus  displacement  is  recorded  to  produce  the  sample 
interferogram.  The  information  in  this  interferogram  enables  one  to  determine  the  sample  thickness,  t. 


MEASUREMENT  METHODOLOGY 

To  understand  how  the  thickness,  t,  of  the  sample  is  derived  from  the  interferogram  obtained 
from  the  reflected-sample  beams,  consider  the  wave  packet  diagram  in  figure  7,  and  the  related 


Figure  7.  Relative  delay  of  front-surface-  and  back-surface- 
reflected  wave  packets  from  sample  by  means  of  Michelson 
interferometer  variable  piezoelectric  transducer  displacement  pt. 
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compuicr-generated  sample  interferogram  in  figure  8.  The  reflected  beams  are  represented  by  wave 
packets  of  width,  lc,  separated  by  the  optical-path-length  difference  2 ngot  between  the  front-surface- 
( fs )  and  back-surface- (bs)  reflected  beams.1  Here,  ng0  is  the  group-refractive  index  of  the  sample 
material,  in  this  case  silicon,  and  has  a  value  equal  to  4.28  for  X0  =  0.83  Hm-2  For  each  of  the  three 
cases  depicted,  two  sets  of  wave  packets  are  shown,  corresponding  to  the  divided  beams  propagating 
in  the  two  arms  of  the  Michelson  interferometer  between  the  beamsplitter  and  the  mirrors. 


TRANSDUCER  DISPLACEMENT 

Figure  8.  Computer-generated  sample  interferogram.  Assumed  values 
areA0  =  0.83  pun,  <5A  =  25  nm,  and  t  =  7.5  pun. 


As  Vp  is  varied  linearly  from  0,  one  set  of  the  reflected  wave  packets  in  the  interferometer  is 
proportionally  delayed  by  an  amount,  2 pt,  with  respect  to  the  other.  At  Vp  =  0,  the  piezoelectric 
transducer  displacement  pt  =  0,  and  ;he  wave  packets  are  coincident  in  space  so  that  they  interfere 
coherently  with  maximum  fringe  visibility.  As  Vp  is  increased,  the  visibility  decreases  as  a  result  of  the 
reduced  overlap  of  the  wave  packets.  (The  source-spectral  width  is  assumed  such  that  the  coherence 
length  is  less  than  the  optical-path  difference  between  the  front-  and  back-surface  wave  packets.)  The 
visibility  reaches  a  minimum  near  pt  =  1/2  ngot  where  there  is  minimum  overlap  between  the  sets  of 
wave  packets.  Then,  the  visibility  increases  as  the  first  surface  reflection  in  one  arm  overlaps  the 
second  surface  reflection  in  the  other  arm.  This  overlap  occurs  at  pt  =  ngot,  where  the  fa-reflected 

’The  optical  path  difference  between  wave  packets  generated  by  a  spectrally  distributed  source  is  proportional  to  the  group 
index  of  refraction,  not  the  phase  index.  For  example,  see  appendix  A  and  Hariharan,  1985. 8 
*This  value  for  ng0  can  be  calculated  using  the  phase  index  values  for  various  wavelengths  found  in  Palik’s  Handbook  of 
Optical  Constant  of  Solids.7 
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beam  in  one  arm  and  the  /y-reflected  beam  in  the  other  arm  are  coherent  because  of  the  compensa¬ 
tion  of  the  optical-path-length  difference  in  the  sample  by  the  optical-path-length  difference  between 
the  arms  of  the  interferometer.  Therefore,  the  difference  in  pt,  measured  between  the  central  maxi¬ 
mum  and  the  side  peak  in  the  interferogram  envelope,  corresponds  to  the  optical  thickness  of  the 
sample,  ngot .  (The  same  result  is  obtained  for  the  opposite  direction  of  piezoelectric-transducer 
displacement,  and  a  maximum  of  the  interferogram  envelope  also  occurs  at  pt  =  -ngot).  Thus,  a 
priori  knowledge  of  the  group  index,  together  with  the  accurate  measurement  of  Pt ,  allows  for  the 
determination  of  film  thickness,  t.  A  more  rigorous  mathematical  derivation  of  the  sample  inter¬ 
ferogram,  including  the  effects  of  dispersion,  is  given  in  appendix  A. 

The  transducer  displacement,  pt,  can  be  determined  by  accurate  calibration  of  the  displacement 
as  a  function  of  Vp,  or  by  the  method  of  fringe  counting  in  the  interferogram.  Fringe-counting  is 
based  on  the  premise  that  each  period  or  fringe  within  the  envelope  of  the  interferogram  corresponds 
to  a  displacement  of  A pt  =  X0/2.  Therefore,  if  sufficient  visibility  exists  over  the  range  of  pt  from  0 
to  ngot,  the  number,  N,  of  fringes  between  the  central  peak  and  one  of  the  side  peaks  can  be 
counted  to  determine  t  by 


~  blXof2  or  t  —  NXq/ 2.ngo  (7) 

For  the  example  shown  in  figure  8,  the  source  was  assumed  to  have  a  symmetric  pure  Gaussian 
spectral  distribution  of  width  6X  =  25  nm,  centered  about  A0  =  830  pm,  and  the  sample  thickness  is 
t  =  7.5  pm.  Counting  the  number  of  fringes  between  the  central  maximum  and  the  secondary  maxi¬ 
mum  yields  N  =  77,  which  from  the  above  relation  gives  t  =  7.47  pm,  in  good  agreement  with  the 
assumed  value.  Note  that  the  uncertainty  in  determining  the  location  of  the  secondary  maximum  is  no 
greater  than  one  fringe,  which  corresponds  to  an  uncertainty  in  t  of  1  IN  or,  in  this  case  1.3  percent. 

Again,  note  that  the  mean  wavelength  of  the  source  distribution  must  be  accurately  known  in 
order  to  calculate  t.  This  can  be  done  by  simultaneously  measuring  the  source-spectral  distribution  via 
the  unused  coupler  port  interfaced  to  a  spectrometer  during  the  recording  of  the  sample  inter¬ 
ferogram. 

In  a  practical  measurement  system,  control  of  the  scanning-piezoelectric  voltage,  spectral 
distribution  measurement,  recording  of  the  interferogram,  determination  of  the  envelope  peaks,  and 
calculation  of  the  sample  thickness,  would  all  be  done  using  an  automated  computer  control  and  data 
analysis  system.  Such  a  system  was  not  used  in  this  work,  but  its  implementation  would  be  straightfor¬ 
ward. 

Fringe-counting  is  an  accurate  means  by  which  the  thickness  of  thin  samples  can  be  measured. 

For  samples  with  thickness  so  large  that  the  optical  path  in  the  sample  is  many  times  greater  than  the 
coherence  length  of  the  source,  the  fringe  visibility  will  be  negligible  in  the  region  between  the  central 
and  secondary  maxima.  Determining  the  occurrence  of  a  fringe  peak  is  not  possible.  In  the  case  of 
large  t  values,  a  supplementary  method  must  be  used.  In  this  regime,  the  accuracy  need  not  be  on  the 
order  of  a  fringe,  since  it  is  only  the  final  etched-sample  thickness  that  must  be  known  with  high 
accuracy.  One  means  to  measure  large  values  of  t,  is  to  mount  one  of  the  Michelson  interferometer 
mirrors  on  a  translation  stage  driven  by  a  stepping  motor.  The  stepping  motor  output  can  easily  be 
controlled  and  monitored  to  provide  patn-length-measurement  accuracies  of  10  pm  over  a  wide  range 
of  translation.  During  the  initial  stages  of  sample  etching,  where  the  thickness  is  on  the  order  of 
1  mm,  the  stepping  motor  would  be  used  to  traverse  the  Michelson  mirror  over  the  distances  neces¬ 
sary  to  cause  the  fs  and  bs  wave  packets  to  overlap.  Then,  in  the  final  stages  of  etching,  where  the 
thickness  is  small  (approximately  10  pm),  the  more  accurate  piezoelectric  transducer  and  the  fringe¬ 
counting  method  would  be  used. 
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EXPERIMENTAL  RESULTS 


Interferograms  were  made  for  various  thin  silicon  CCD  samples  with  a  nominal  thickness  of  8  gm. 
The  ELED  temperature  and  current  were  held  constant  at  25  °C  and  10  mA.  The  source  spectral 
distribution  under  these  conditions  is  shown  in  figure  9.  The  mean  wavelength  is  X0  =  0.847  pm,  and 
the  spectral  width  is  approximately  6X  =*  20  nm  or  do  —  0.028  pm"1 .  Some  periodic  structure  and 
asymmetry  are  also  present,  but  these  effects  proved  to  be  inconsequential  as  discussed  in  the  follow¬ 
ing  section. 


WAVELENGTH  (nm) 


Figure  9.  Spectral  distribution  for  ELED  used  in  the  experiments  operating  at  10  mA. 
Mean  wavelength  is  847  nm  and  spectral  width  is  approximately  20  nm. 


The  output  beam  of  the  fiber  coupler  was  focused  onto  a  small  region,  about  100  pm  in  diameter, 
at  the  center  of  the  sample.  The  center  was  found  by  displacing  the  sample  along  orthogonal  axes  in  a 
plane  normal  to  the  beam  and  locating  the  positions  of  the  sample  edges  by  monitoring  the  intensity 
of  the  reflected  beam. 

A  typical  interferogram,  obtained  by  photographing  the  oscilloscope  trace  of  the  detector-amplifier 
voltage,  Vd,  versus  the  transducer-drive  voltage,  Vp ,  is  shown  in  figure  10.  In  this  photograph,  the 
number  of  fringes  can  be  counted  visually  to  yield  N  =  79,  which,  using  equation  (7),  corresponds  to  a 
sample  thickness  of  7.8  pm  at  the  center.  Numerous  silicon  CCD  samples  were  measured  and  the 
results  are  tabulated  in  table  1. 
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An  attempt  was  made  to  determine  if  any  variation  of  sample  thickness  exists  as  a  function  of 
position  on  the  sample.  Again,  the  edges  of  the  2  mm  x  2  mm  sample  were  located  and  a  matrix  of 
points,  equally  spaced  at  500  )im  intervals  within  this  area,  was  probed  by  the  output  beam  of  the 
coupler.  The  results  for  one  particular  sample  are  shown  in  table  2,  which  shows  that  the  maximum 
variation  of  thickness  over  the  sample  is  0.3  jj.m. 


TRANSDUCER  DRIVE  VOLTAGE.  Vp 

Figure  10.  Actual  sample  interferogram  obtained  by 
photographing  oscilloscope  recording  of  detector-amplifier- 
output  voltage  versus  piezoelectric-transducer-ramp  voltage. 
Number  of  fringes  from  central  maximum  to  secondary 
maximum  is  N  =  79,  corresponding  to  a  sample  thickness 
of  t  =  7.8  jim. 


Table  1.  Results  of  sample  thickness  measurement. 


Hughes  CCD  Sample  ft 

N 

t  (M-m) 

7153AA20R1-6F 

87 

8.6 

7154AA19R1-6J 

78 

7.7 

7153AA20R1-5J 

86 

8.5 

7154AA5R1-4D 

77 

7.6 

7154AA3R1-5E 

84 

8.3 

7154AA5R1-5A 

80 

7.9 

7153AA09R2-7C 

76 

7.5 

7153AA09R2-4G 

81 

8.0 

7153AA14R1-2C 

77 

7.6 

7153AA14R1-2J 

79 

7.8 

7154AA3R1-6A 

84 

8.3 
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Table  2.  Sample  thickness  variation  with  position. 


x  (pm) 

-500 

0 

+500 

+500 

7.7 

7.9 

7.7 

y  (Mm) 

0 

7.6 

7.8 

7.8 

-500 

7.7 

7.6 

7.6 

DISCUSSION  AND 
CONCLUSIONS 

A  method  has  been  described  for  the  precision  measurement  of  thin  silicon  CCD  samples  using  a 
fiber-optic-coupler  probe,  an  edge-emitting  LED  (ELED),  and  a  Michelson  interferometer.  The  meas¬ 
urement  precision  is  one  optical  fringe,  or  about  0.1  pm.  The  fiber-optic  probe  can  do  remote  in  situ 
measurements  and  variable  position,  localized  thickness  measurements.  The  flexibility  of  the  fiber,  and 
the  small  size  and  chemical  inertness  of  the  fiber/GRIN  rod  lens  assembly,  enable  the  probe  end  to 
be  configured  in  a  probe  housing  that  could  be  placed  close  to  the  sample  within  the  etching  chamber 
while  the  sample  is  processed.  In  this  way,  continuous  sample-thickness  measurements  could  be  made, 
and  the  data  used  to  control  the  etching  parameters. 

While  the  basic  operation  of  the  equipment  and  measurement  methods  have  been  successfully 
demonstrated,  several  limitations  and  potential  future  modifications  need  to  be  considered.  As  men¬ 
tioned  previously,  the  spectral  distribution  of  the  source  determines  several  basic  features  of  the 
sample  interferogram.  The  source  width  (hence,  the  coherence  length)  must  be  such  that  the  fringe 
visibility  goes  through  a  distinct  minimum,  while  maintaining  adequate  visibility  for  fringe-counting,  as 
the  Michelson  interferometer  is  scanned  to  generate  the  interferogram.  For  samples  with  thickness 
values  on  the  order  of  5  to  10  pm,  the  source  must  be  either  an  ELED  or  a  superluminescent  diode 
(SLD).  A  suitable  SLD  was  not  available  for  this  work,  and  therefore  the  system  could  not  be  evalu¬ 
ated  with  this  type  of  source.  (The  SLD  offers  potentially  higher  optical  powers  than  the  ELED,  which 
could  be  important  in  certain  applications,  such  as  in  situ  measurements,  where  the  optical  power  of 
the  reflected  sample  beam  might  be  small.) 

Either  type  of  source  can  exhibit  both  spectral  structure  and  asymmetry.  It  is  important  that  the 
structure  be  minimized  so  that  the  source  interferogram  components  do  not  interfere  with  the  sample 
interferogram  secondary  maximum.  In  this  work,  the  ELED  source  spectral  properties  were  carefully 
investigated.  Source  spectra  measured  at  a  constant  temperature  of  25 °C  for  various  drive  currents  are 
shown  in  figure  11.  Note  that  there  is  significant  source  structure  and  asymmetry  for  the  higher  oper¬ 
ating  currents,  and  that  the  mean  wavelength  shifts  to  higher  values  as  the  current  is  decreased.  To 
minimize  the  effects  of  source  spectrum  periodic  structure  and  asymmetry  discussed  in  the  previous 
sections,  the  ELED  was  operated  at  a  current  of  10  mA.  (Even  at  this  low-operating  current  there 
was  adequate  light  energy  to  provide  a  large  signal-to-noise  ratio  (SNR)  for  the  system  detector  out¬ 
put.)  Figure  9  is  an  expanded  plot  of  the  spectral  distribution  for  an  operating  current  of  10  mA.  It 
can  be  seen  that  even  at  this  low-operating  current  some  structure  is  present  with  a  period  of 
AX  =  5.4  nm  (Ao  =  0.0075  pm"1). 
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Figure  11.  ELED  source  spectra  at  25°C  for  various  operating  currents.  Presence  of 
structure  and  asymmetry  increase  as  current  increases,  while  mean  wavelength  shifts  to 
smaller  values. 

The  effect  of  this  structure  is  minimal  in  the  measurement  of  thin  silicon  films.  First,  the  period  is 
such  that  a  contribution  to  the  source  interferogram  would  be  made  at  a  path  difference  of 
A p  =  1/A  a  =  133  jim,  which  is  much  greater  than  the  sample  round-trip  optical-path  length  of 
2 ngot  —  60  jim.  Second,  the  amplitude  of  the  periodic  structure  is  much  less  than  the  Gaussian 
component  amplitude  and,  therefore,  its  contribution  to  the  interferogram  will  be  correspondingly 
small.  These  conclusions  are  confirmed  by  the  10-mA  source  interferogram  shown  in  figure  12  for  a 
range  of  A p  from  -10  jim  to  110  jim. 

The  mean  wavelength  must  also  be  considered  in  the  selection  of  an  appropriate  optical  source 
and  depends  on  its  material  composition.  Both  GaAlAs  and  InGaAsP  ELEDs,  operating  in  the  spec¬ 
tral  regions  near  0.83  jim  and  1.3  jim,  respectfully,  are  commercially  available. 

Each  fringe  in  the  interferogram  corresponds  to  an  optical-path  difference  of  A0/2,  and,  therefore, 
the  shorter  the  wavelength,  the  greater  the  precision  in  the  measurement  of  the  sample  thickness.  In 
this  respect,  the  GaAlAs  source  with  X0  =*  0.83  jim  would  be  better  than  an  InGaAsP  source.  On  the 
other  hand,  the  absorption  of  light  is  greater  in  silicon  for  shorter  wavelengths,  and,  therefore,  for 
relatively  thick  samples,  the  back-surface-reflected  beam  would  be  too  strongly  attenuated  to  contrib¬ 
ute  to  the  interferogram. 

The  value  of  the  interferogram  secondary  maximum,  relative  to  the  maximum  at  0  path-length 
difference,  is  related  to  the  ratio  of  the  fields  resulting  from  the  back-  and  front-surface  reflections.  A 
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ratio  of  10  percent  is,  somewhat  arbitrarily,  chosen  as  the  minimum  value  that  will  result  in  a  clearly 
discernible  secondary  maximum  for  fringe-counting. 
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TRANSDUCER  DRIVE  VOLTAGE.  Vp 

Figure  12.  Source  interferogram  for  ELED  operating  at 
10  mA.  Minimal  periodic  structure  is  indicated  by  lack 
of  significant  interferogram  amplitude  other  than 
central  envelope. 


The  back-surface  field  is  reduced  with  respect  to  the  first  surface  field  by  three  factors  in  addition 
to  the  sample  absorption.  First,  the  coupling  of  the  reflected  beams  back  into  the  fiber  will  be  differ¬ 
ent  due  to  the  unequal  distances  from  the  principal  plane  of  the  GRIN  rod  lens.  (In  practice,  the 
relative  position  of  the  lens  and  the  sample  is  adjusted  so  that  the  height  of  the  secondary  maximum 
is  maximized  for  a  given  film  thickness.)  Second,  the  back-surface  reflection  is  transmitted  twice 
through  the  sample/air-first  surface,  and,  therefore,  the  back-surface-field  amplitude  is  reduced  by  a 
factor  of  t2  relative  to  the  front-surface  field.  Here,  t  is  the  coefficient  of  the  electric  field  transmis¬ 
sion,  and  for  silicon  in  t2  =*  0.68.  Third,  while  the  front  surface  interface  is  with  air,  the  back  surface 
has  on  it  the  CCD  circuitry,  which  consists  of  thin  layers  of  polymers,  metals,  and  oxides.  These 
materials  are  not  necessarily  uniform,  nor  highly  reflecting,  with  the  result  that  the  back-surface 
reflection  is  less  than  the  front-surface  silicon/air  reflection.  The  reduction  of  the  back-surface  reflec¬ 
tion  due  to  these  factors,  and  the  absorption  by  the  material  itself,  determine  the  source  wavelength 
limitations  on  the  sample  thickness  that  can  be  measured. 

If  the  coupling  efficiency  back  into  the  fiber  is  assumed  to  be  equal  for  the  front-  and  back- 
surface-sample-reflected  beams,  that  the  reflection  coefficients  are  equal  for  the  two  surfaces,  and 
taking  into  account  the  reduction  in  the  back-surface  field  by  the  interface  transmittance  factor,  t2, 
then  a  limit  on  the  maximum  thickness  of  the  sample  can  be  determined  due  to  absorption.  The  ratio 
of  the  interferogram  envelope  secondary  and  primary  maxima  is 


EfsEbs/(E}s  +  Els) 


(8) 


14 


where 


Ebs/Efs  =  t2exp(-  at)  (9) 

and  a  =  is  the  absorption  coefficient  of  silicon.  For  X  =  0.83  pm,  a  =  700  cm-\  while  for 
\  =  1.3  pm  silicon  is  transparent  (a  -*  0).8  Equations  (8)  and  (9)  can  be  used  to  calculate  the  maxi¬ 
mum  sample  thickness  that  would  yield  a  secondary  interferogram  envelope  maximum  equal  to 
10  percent  of  the  primary  maximum.  The  result  is  *max-  27  pm  for  X  =  0.83  pm. 

Note  that  in  figure  10,  the  ratio  of  the  interferogram  maxima  for  a  7.6  pm  sample  is  approxi¬ 
mately  10  percent.  This  is  an  indication  that  the  assumptions  regarding  beam-to-fiber  coupling  efficien¬ 
cies  and  equal-interface-reflection  coefficients  are  not  necessarily  valid.  In  practice,  for  a  sample 
thickness  greater  than  10  pm,  the  ratio  will  always  be  less  than  10  percent  at  0.83  pm  due  to 
increased  absorption;  therefore,  the  longer  wavelength  InGaAsP  source  would  be  more  suitable. 
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GLOSSARY  OF  SYMBOLS  AND  ABBREVIATIONS 


bs 

CCD 

E(o) 

Es(o) 

Efs 

Ebs 

ELED 

fs 

GRIN 

HeNe 

lc 

LED 

mA 

n 

n8 

ngo 

no 

"o' 

N 

nm 

p.m 

Py 

P2 

Pt 

SLD 

SNR 


back  surface 

charge-coupled  device 
interferometer-output-optical  field 

source-optical  field 

optical  field  from  sample-front  surface 
optical  field  from  sample-back  surface 
edge-emitting,  light-emitting  diode 
front  surface 
graded  index 
Helium  Neon  (laser) 
source-coherence  length 
light-emitting  diode 
milliampere 

phase  index  of  refraction  of  the  sample  material 
group  index  of  refraction  of  the  sample  material 
group  index  of  refraction  of  the  sample  material  at  a  =  o„ 
phase  index  of  sample  material  for  a  =  aQ 
derivative  of  phase  index  for  a  =  a0 

number  of  interferogram  fringes  between  primary  and  secondary  maxima 

nanometer 

micrometer 

optical  path  length  of  (first)  wave  packet 
optical  path  length  of  (second)  wave  packet 
differential  length  between  the  interferometer  arms 
superluminiscent  diode 
signal-to-noise  ratio 
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t  thickness  of  the  sample  to  be  measured 

t  optical  field  transmittance  at  sample/air  interface 

Vd  photodetector  output  voltage 

Vp  voltage  applied  to  interferometer  transducer 

a  absorption  coefficient  of  silicon 

<5ct  width  of  spectral  energy  distribution  in  wavenumbers 

<5A  spectral  width  of  source  in  wavelength  units  (nm  or  p.m) 

A p  path-length  difference  between  optical  wave  packets 

AA  period  of  spectral  structure  in  wavelength  units 

Act  period  of  spectral  structure  in  wavenumber  units 

A  wavelength  of  optical  field  (nm  or  (im) 

X0  mean  wavelength  of  source  distribution 

ct  wavenumber  of  optical  field  (p.m_1) 

°o  mean  wavenumber  of  source  distribution 


17 


APPENDIX  A 


THEORETICAL  DERIVATION  OF  THE  INTERFEROGRAM  FOR  A 
SYMMETRIC  MICHELSON  INTERFEROMETER  WITH 
A  TWO-BEAM  INPUT  FROM  A  DISPERSIVE  MEDIUM 

In  this  appendix,  a  more  detailed  derivation  of  the  sample  interferogram  is  presented.  Again  the 
interferometer  is  a  symmetric  Michelson  interferometer  with  a  two-beam  input  from  the  sample;  but 
now,  the  effects  of  material  dispersion  are  included.  As  in  the  previous  treatment  of  the  source 
interferogram,  the  inputs  to  the  interferometer  are  assumed  to  have  the  same  polarization,  can  be 
treated  as  emanating  from  point  sources,  and  the  beams  have  no  angular  misalignment.  The  use  of 
the  fiber  coupler  input  assures  that  all  of  these  assumptions  are  valid. 

It  will  be  shown  that  for  a  typical  sample  thickness  of  the  order  of  10  pm,  the  effects  of  material 
dispersion  can  be  neglected,  and  that  the  previous  results  regarding  the  features  of  the  sample  inter¬ 
ferogram,  which  were  based  on  the  less  qualitative  wave-packet  overlap  treatment,  are  completely 
valid. 

Referring  to  figure  6,  the  continuous  optical  output  field  from  the  source,  represented  by  Es(o),  is 
incident  on  the  sample  via  the  bidirectional  coupler.  This  beam  is  partially  reflected  from  the  front 
and  back  surfaces  with  a  relative  phase  difference  between  the  reflected  fields  of  Ajiant ,  correspond¬ 
ing  to  the  optical-path  length  in  the  material.  Here  n  is  the  phase  index  of  the  material  evaluated  at 
the  particular  wavenumber  o .  Due  to  the  partial  transmittance  of  the  optical  field  at  the  sample 
material  boundaries,  the  amplitude  of  the  back-surface-reflected  beam  relative  to  the  front-surface 
reflection  is  t2,  where  t  is  the  transmittance  of  the  sample/air  interface.  The  front-  and  back-surface- 
reflected  beams  are  coupled  into  the  Michelson  interferometer  via  the  GRIN  rod  lens  and  the  fiber 
coupler.  The  differential  optical-path  length  between  the  two  arms  of  the  interferometer  is  2 pt,  and  is 
varied  by  means  of  the  piezoelectric  transducer. 

Using  the  complex  field  notation,  and  omitting  the  explicit  time  dependence  of  the  fields,  the 
output  field  from  the  interferometer  is  expressed  as 

E{a)  =Re{U{o)}  (A-l) 

where  U[o)  is  the  complex  field  amplitude  as  a  function  of  a, 

U(o)  =  Es(o)  [  1  +exp(/4;ropO][l  +  t2exp0'4jrom)] 

=  £i(a)[l  +  t2expO'4jranO  +  exp(/4jropf)  + 12  exp(/47rop(  +  nr)  ]  (A-2) 

=  UX  +  Ut  +  Up+U2 


(Note  that  (/]  =Es(o),  and  that  common  phase  and  amplitude  factors  have  been  omitted.) 

In  this  expression,  it  is  understood  that  n  is,  in  general,  a  function  of  a .  Also,  it  has  been 
assumed  that  the  mirror  reflectances  are  both  equal  to  1,  and  that  the  beamsplitter  is  lossless.  The 
time-averaged  intensity  for  a  wavenumber  component  a  is  then  expressed  as 


A-l 


(A-3) 


1(a)  =  \U(o)\2  =  UU  * 

=  |C/1|2+K/,|2  +  |t^|2  +  |l/2|2  +  (l/1l/t*  +  (/1  •  Ut)  +  (UPU2  *  +  Up  *  U2) 

+  (UxUp  *  +  (/,*  Up)  +  (t/,l/2  ’  +  U,  *  t/2) 

+  (t/,t/2  *  +  U,  *  t/2)  +  *  +  ut*up 

In  the  above  equation,  the  terms  have  been  grouped  according  to  the  constant  terms  on  the  first  line, 
which  are  not  a  function  of  pt,  and  the  remainder  of  the  terms,  which  are  variable.  This  expression 
can  be  further  simplified  by  using  the  relation  (Z+Z*)  =  2Re{Z),  where  Z  is  any  complex  number,  and 
by  noting  that 

|(/1|2  =  |l/p|2  =  |£i(a)|2 

\U'\2  =  \U2\2=l2\Es(o)\2 
(UpU2  *  +  Up  *  U2)  =  UXU,  *  +  Ur*  U, 
and  (UtU2  •  +  U,*U2)=  1 2(U2UP  *  +  U,  *  Up) 


Using  these  relations, 

/(a)  =  2 [ ( 1  +t2)|£'^(a)|2  +  2  Re {(/,(/,}  +  (1  + 12)  Re{UxUp)  +  /?e{U1t/2}+  Re{UtUp*]  (A-4) 

The  first  three  terms  are  constants  with  respect  to  pt,  and  contribute  to  the  constant  intensity 
present  in  the  output.  The  last  three  terms  contain  the  information  related  to  the  sample  thickness  t. 

The  total  intensity  I(pt)  of  the  output  beam  is  obtained  by  integrating  the  function  1(a)  over  the 
entire  spectral  range  of  the  source, 


I  (a)  da 


Denoting  the  integrals  of  the  three  constant  terms  in  the  expression  for  I  (a)  as  /(“>),  the  expres¬ 
sion  for  I(pt)  becomes 


f(p<)  =/(»)+  2  Re 


J'gUi  U2do  + 


(A-5) 


Because  of  the  limited  spectral  range  of  the  source,  the  limits  of  integration  have  been  chosen  to  be 
±  ®  for  mathematical  convenience.  The  notation  for  /(<»)  follows  by  recognizing  that  as  pt  -*  °°,  the 
three  nonconstant  terms  -*  0,  and  therefore,  the  first  three  terms  simply  represent  the  value  of  I(pt) 
at  pt  -*  oo. 

At  this  point,  a  specific  functional  form  for  the  source  energy  spectral  distribution,  |£i(a)|2,  is 
chosen  which  represents  typical  high-radiance-limited-spectral  width  sources,  such  as  light-emitting 
diodes  (LEDs)  and  superluminescent  diodes  (SLDs).  These  sources  can  be  represented  by  a  normal¬ 
ized  Gaussian  spectral  distribution  of  the  form 


|£<(o)|2  =  (l/<5cr)  exp  {- n(o - Oo)1  / (do)2}  . 


(A-6) 


A-2 


To  evaluate  the  integrals  in  equation  (A-5),  a  change  of  variables  is  made  so  that  the  integrals  are 
symmetric  about  o0.  Letting  s  =  a-a0,  and  setting  the  limits  of  integration  on  s  to  be  ±  00  (this  is 
valid  since  o0»5a),  each  of  the  three  integrals  in  equation  (A-5)  can  be  put  in  a  standard  form  that 
can  be  found  in  integral  tables.  The  general  relation  from  reference  1  is 


+  00 


{- [as2 +  2bs +  c]ds  =  (n/a)*  exp  {(b2-ac)/a}  . 


(A-7) 


Here,  a  and  b  are  the  appropriate  constants  for  each  integral.  Using  equations  (A-5),  (A-6),  and 
(A-7),  the  total  intensity  is  then 

Hpt)  -■  /(°o)  +  2  Re(  (1  + t2)  exp[;4jror0/>*]  exp[- ;r(p()2(2<5a)2j 


2t2  ex\>[iAjio0(pt  +  n0t)  ] 


[  1  -  in„'/(2<5a)2]2 


i 

2lT 


exp  {-  ji(pt  +  ngot)2(26o)2/[l  -  in0't{2do)2\} 


(A-8) 


2t2  exp[-  iAnooipt  -  n0t ] 


[1  -  in0't{25o)2 


-p^-exp  { -  ji(pt  -  ng0t)2(26o)2/[\  -  in0't(2da)2]}^ 


To  obtain  this  result,  it  has  been  assumed  that  the  phase  index  of  refraction,  which  is  a  function  of 
the  wavenumber,  can  be  expanded  in  a  Taylor  series  near  o0, 


n=n(o)  =  n(a0)  +  (o-o0)(dn/do)\o0 


=  n0+  ( o-o0)n0 '  . 


Also,  the  notation  ngo  has  been  used  for  the  value  of  the  group  index  at  a  =  o0. 
ngo  =  ^(o)|ao  =  [n(o)  +o(dn/da))|ao  . 

Using  equation  (A-8),  the  features  of  the  interferogram,  which  were  derived  previously  using  the 
wave  packet  overlap  analysis,  are  confirmed.  First,  the  interferogram  function  is  an  oscillatory  function 
with  an  envelope  that  consists  of  three  Gaussian  functions  centered  at  pt  =  0  and  ±  ngot.  Also,  the 
shape  of  the  envelope  function,  and  the  phase  of  the  oscillation,  depend  on  the  source  mean 
wavenumber  and  spectral  width,  the  sample  thickness,  and  the  phase  and  group  refractive  indices  of 
the  sample. 

To  this  point,  the  derivation  of  the  interferogram  function  is  valid  for  a  dispersive  sample  medium 
of  any  thickness.  The  contribution  of  dispersion  appears  whenever  the  factor  n0'  is  present,  and  the 
magnitude  of  the  effect  increases  with  increasing  sample  thickness.  It  can  be  seen  that  dispersion 
affects  both  the  phase  and  amplitude  of  the  interferogram  function.  The  amplitude  reduction  of  the 
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secondary  maxima  is  not  significant  as  long  as  there  is  sufficient  signal  to  allow  for  peak  height  deter¬ 
mination  for  fringe-counting.  The  phase  is  affected  by  both  the  addition  of  constant  phase  terms,  and 
by  a  chirping,  or  nonlinear  variation  of  phase  with  p/. 

To  assess  the  effect  of  dispersion,  it  is  helpful  to  consider  the  magnitude  of  these  factors  for  a 
particular  experimental  situation  corresponding  to  a  practical  sample  measurement.  The  source  used  in 
these  experiments  was  chosen  to  be  an  ELED  with  a  mean  wavelength  of  approximately  0.85  pm  and 
a  spectral  width  of  approximately  20  nm  (see  figure  9).  These  values  correspond  to  values  of 
o0  =  1.18pm"1  and  6a  =  0.028  pm'1  .  The  sample  is  a  silicon  CCD  with  a  thickness  on  the  order  of 
10  pm.  The  values  of  the  group  index  and  the  derivative  of  the  phase  index  at  Pt  can  be  derived 
from  published  data  in  reference  7,  and  are  found  to  be  ngo  =  4.28  and  n0'  =  0.5  pm. 

For  these  values  of  pt,  6a,  t  and  n0'  it  is  found  that  the  contributions  to  the  interferogram  due  to 
material  dispersion  can  be  neglected.  For  example,  the  constant  phase  factor  is 

<p  =  tan"1(2n0'r((5(j)2)  =  tan"‘(0.008)  =  0.008  rad 

Similar  negligible  values  are  obtained  for  the  amplitude  reduction  factor  and  the  chirp  phase  term. 

It  should  be  noted  that  these  negligible  values  are  for  small  values  of  t,  and  that  for  larger  values, 
these  effects  should  be  included.  For  example,  if  the  maximum  tolerable  phase  offset  due  to  material 
dispersion  is  chosen  to  be  10  percent  for  either  the  second  or  third  terms  in  the  interferogram,  then 
the  corresponding  maximum  sample  thickness  is  =  50  pm. 

If  the  sample  thickness  is  small  enough  that  the  material  dispersion  effects  can  be  ignored,  then 
the  interferogram  function  can  be  simplified  to  give  (omitting  the  t  dependent  amplitude  factors) 

I(pt)  =/(«)  +cos  [4jio0pt}  expf-Tr^f)2^)2] 

+  cos  (4w0(p,  +  n0/))  exp[-^(p(  +  ni?o02(25a)2]  (A-9) 

+  cos  {47io0{pt  -n0t)}  exp[-7r(p/  -  ngot)2(26a)2] 

Equation  (A-9)  is  the  mathematical  representation  of  the  result  that  was  derived  qualitatively  using  the 
wave  packet  approach,  and  which  is  represented  by  the  interferogram  shown  in  figure  8. 
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